Ballistic heat transport of quantum spin excitations as seen in SrCu02 
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Fundamental conservation laws predict ballistic, i.e., dissipationless transport behaviour in one- 
dimensional quantum magnets. Experimental evidence, however, for such anomalous transport has 
been lacking ever since. Here we provide experimental evidence for ballistic heat transport in a 
S = 1/2 Heisenberg chain. In particular, we investigate high purity samples of the chain cuprate 
SrCu02 and observe a huge magnetic heat conductivity Kmag. An extremely large spinon mean 
free path of more than a micrometer demonstrates that Kmag is only limited by extrinsic scattering 
processes which is a clear signature of ballistic transport in the underlying spin model. 

PACS numbers: 75.40.Gb, 66.70.-f, 68.65.-k, 75.10.Pq 



The integrability of the one-dimensional (ID) antifer- 
romagnetic S = 1/2 Heisenberg chain implies highly 
anomalous transport properties, in particular, a divergent 
magnetic heat conductivity Krnug at all finite tempera- 
tures T.—"— This truly ballistic heat transport suggests 
anomalously large life times and mean free paths of the 
quantum spin excitations and renders ID quantum mag- 
nets intriguing candidates for spin transport and quan- 
tum information processingi^^— However, despite the rig- 
orous prediction, experimental evidence for ballistic heat 
transport in quantum magnets is lacking. Nevertheless, 
promising large Kmag has been observed in a number of 
cuprate compounds which realize ID S = 1/2 Heisen- 
berg antiferromagnets^rJ^ with the spin chain material 
SrCu02 being a prominent exampl o^^i^^ although a quan- 
titative analysis of Kmag has always been difficult there 
since the phononic and magnetic heat conductivities are 
of similar magnitude at low temperature. Such exper- 
imental Kmag is always finite since extrinsic scattering 
processes due to defects and phonons are inherent to all 
materials and mask the intrinsic behavior of the chain. 
Formally it seems reasonable to account for the extrinsic 
scattering via a finite Kmag ^ DthT, where r is a relax- 
ation time, and Dth represents the thermal Drude weight 
which (multiplied by a delta function at zero frequency) 
describes the intrinsic heat conductivity. In fact, it was 
thereby possible to identify the expected low-T linear- 
ity of Dth (7^) in the case of a "dirty" spin chain material 
where a high density of chain defects generate a large 
T- independent scattering rate 1/t^ 

In this paper we examine the heat conductivity of 
SrCu02 which is considered an excellent realization of 
the 5=1/2 Heisenberg chainji^Ti^ Our samples of ex- 
traordinary purity allow an unambiguous separation of 
the phononic and magnetic contributions to the thermal 
conductivity. This yields the by far highest Kmag ob- 
servediSiii^ until now. Our analysis reveals a remarkable 
lower bound for the low-T limit of the mean free path 
I 



increasing temperatures Kmag is increasingly supressed 
due to spinon-phonon scattering which is the dominant 
extrinsic scattering mechanism in this material. 

We have grown large single crystals of SrCu02 by 
the traveling solvent floating zone method,— where the 
feed rods were prepared using the primary chemicals 
CuO and SrCOa with both 2N (99%) and 4N (99.99%) 
purity. Cuboidal samples with typical dimensions of 
(3 X 0.5 X 0.5) mm-^ were cut from the crystals, with the 
longest dimension parallel to the principal axes. Four- 
probe measurements of the thermal conductivity k were 
performed in the 7-300K range^^ with the thermal cur- 
rent along the a, b, and c-axes (k^, k;,, and Kc respec- 
tively) for both the 2N and the 4N samples. 

The main structural element in SrCu02 is formed 
by Cu02 zig-zag ribbons, which run along the crys- 
tallograpliic c-axis (see inset Fig. [ij. Each rib- 
bon can be viewed as made of two parallel corner- 
sharing Cu02 chains, where the straight Cu-O-Cu 
bonds of each double-chain structure result in a 
very large antiferromagnetic intrachain exchange cou- 
pling J/kB « 2100 - 2600 K of the 5 = 1/2 
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of more than a micrometer. Thus our data pro- 
vide striking evidence that the intrinsic heat transport 
of the S=l/2 Heisenberg chain is indeed ballistic. With 



weaker interchain coupling 

and presumably quantum fluctuations prevent three- 
dimensional long range magnetic order of the system at 
T > Tat « 1.5-2 K « lO^V/fcs KM^ Hence, at signif- 
icantly higher T the two chains within one double chain 
structure can be regarded as magnetically independent. 
In fact, low-T (12 K) inelastic neutron scattering spectra 
of the magnetic excitations can be very well described 
within the 5=1/2 Heisenberg antiferromagnetic chain 
mo del 1^ 

Fig. [T] presents our results for Ka and Kc of SrCu02 for 
both 2N and 4N purity as a function of T. We first de- 
scribe the data for 2N purity which are in good agreement 
with earlier results by Sologubenko et alJ^ A pronounced 
low-T peak at -18 K with K^ax « 215Wm"iK"i is 
found for 2Nj i-e., perpendicular to the chains. This 
peak and a — T~^-decrease at T > 150 K towards a 
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Figure 1: (Color online) Ka and Kc of SrCu02 for different 
purity values. Closed (open) symbols represent c-axis (a-axis) 
data, circles (diamonds) correspond to 4N (2N) purity. Inset: 
crystal structure of SrCu02. The symmetry is Cmcm with 
lattice constants a = 3.56 A, 6 = 16.32 A, c = 3.92 A.^^ 



small value at room temperature (^ 6 Wm^^K^^) repre- 
sent the characteristic T-dependence of phonon-only heat 
conductivity Kph- The peak originates from two compet- 
ing effectsi^l at low T, a weakly T-dependent phonon 
mean free path /ph and a rapidly increasing number of 
phonons cause Kph to increase strongly. At higher T, 
the exponentially rising number of phonon-phonon umk- 
lapp processes increasingly shortens Zph, which causes the 
decrease of Kp^. A similar low-T peak (at ^ 20 K) 
is also present in Kc,2N (parallel to the chains). It is 
however larger (^max ~ 335 Wm~-'^K~^) and exhibits a 
distinct shoulder at the high-T edge (T > 40 K) of 
the peak. Kc.2N decreases at higher T, but remains 
much larger than Ka.2N and even at room temperature 
Kc.2N ^ 40Wm^^K^^. The apparent large anisotropy, 
together with the unusual T-dependence of Kc,2N, is the 
signature of a large magnetic fraction of Kc,2N over a large 
r-range3i^ 

We now turn to the new data which have been ob- 
tained for the high-purity compound. The heat trans- 
port perpendicular to the chains (Ka,4N) is slightly en- 
hanced as compared to Ka.2N (^max ~ 235Wm~^K~^) 
which reflects a somewhat reduced phonon-defect scat- 
tering. However, a much more drastic and unexpected 
large effect of the enhanced purity is observed in the 
heat transport parallel to the chains, Kc.4n- Instead of a 
narrow low-T peak and a shoulder as observed in Kc.2N, 
a huge and broad peak centered at ~28 K is present 
in k;c,4N (^max ~ 830Wm~^K~^) which exceeds Kc,2N 
at T < 70 K by more than a factor of 2. Also at 
70K < T < 300 K we observe Kc,4n > Kc.2n, where 



interestingly both curves approach each other and at 
T > 200 K exhibit almost the same T-dependence. 

Without further analysis some clear-cut conclusions 
can be drawn. First, the extraordinary enhancement 
of Kc upon the improvement of the material's purity in 
contrast to a concomitantly negligibly small one in Ha, 
straightforwardly implies that the enhancement primar- 
ily concerns the magnetic heat conductivity Kmag which 
is present in Kc only. Second, the extreme low-T sensi- 
tivity to impurities of Kmag suggests that spinon-defect 
scattering is the dominating process which relaxes the 
heat current in this regime. Third, upon rising T, the 
spinon-defect scattering is increasingly masked by a fur- 
ther scattering process which leads to Kc,2N and Kc,4N 
being very similar at T > 200 K. The most reasonable 
candidate for this process is spinon-phonon scattering, 
since the only thinkable alternative, i.e. spinon-spinon 
scattering, is negligible^"— in this T-regime. 

A further analysis of the data requires a reliable sep- 
aration of the total measured k into all relevant contri- 
butions which normally add up. Since electronic con- 
tributions can be excluded in this electrically insulating 
material, it seems natural to assume that the measured 
Kc is just the sum of Kmag and a phononic background 
i^ph,cr'— where the latter can be approximated by the 
purely phononic heat conductivity perpendicular to the 
chains, Ka ~ ki, (see inset of Fig. [5]). The thus obtained 
'^mag = Kc — Ka for the 2N and the 4N samples are shown 
in Fig. [51 At T < 35 K, i.e., in the vicinity of the peak of 
'tph.c, errors become large and we disregard the data in 
this range for further analysis. For higher T we account 
for a possible uncertainty of ±30% in Kph.c- Note that, 
in the case of the 4N compound, possible errors in Kmag 
are rendered small because obviously Kmag,4N S> HaAisi- 

^mag of the 4N sample exhibits a sharp peak at 
^37 K with an extraordinary maximum value of about 
660Wm"^K"^, which is more than a factor of 3 higher 
than the largest reported Kj^ag'^'^ The peak is followed 
by a strong decrease upon raising T. Similar to the afore 
described typical T-dependence of a clean phononic heat 
conductor, the overall T-dependence of Kmag suggests 
that, in a simple picture, two competing effects deter- 
mine Kmag- The low-T increase of Kmag is consistent 
with a regime where the effect of scattering processes is 
weakly T-dependent since Dth is expected to increase lin- 
early with Ti^"— The strong decrease at higher T is then 
the result of the increasing importance of spinon-phonon 
scattering. Kmag of the 2N sample is qualitatively very 
similar. However, the absolute value at the peak is much 
lower (^ 172 W/mK) and the peak's position is shifted to 
a higher T (~ 55 K). Similarly to the Kc data, at higher 
T, the 4N curve approaches that of the 2N sample. The 
latter is consistent with the earlier notion that spinon- 
phonon scattering is dominant at high T, while the dif- 
ferences at low T suggest that spinon-phonon scattering 
freezes out, upon decreasing T, rendering spinon-defect 
scattering increasingly important. 

We analyze Kmag quantitatively by extracting the 
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Figure 2: (Color online) Kmag of SrCu02 for different purities. 
Open symbols represent low-T Kmag which is disregarded in 
the further analysis. The shaded areas show the uncertainty 
of the estimation of Kmag due to the phononic background. 
Inset: Ka and ki, perpendicular to the chain for both purities. 



spinon mean free path Imug according to^ — 



where Ns = 4/a6 is the number of spin chains per unit 
area. As can be inferred from Fig. [31 /mag of both samples 
show a strong decrease with increasing T, which directly 
reflects spinon-phonon scattering becoming increasingly 
important. Both curves are very similar, but clear differ- 
ences are present at low T, where /mag of the 2N sample 
is somewhat lower, in accordance with a higher spinon- 
defect scattering. We evoke Matthiesen's rule to model 
the T-dependence of Zmag and to account for both scatter- 
ing processes viz. /"^g = lo^^^sp- Here, Iq describes the 
T- independent spinon-defect scattering whereas lsp(T) 
accounts for the T-dependent spinon-phonon scattering. 
For the latter, we assume a general umklapp process with 
a characteristic energy scale ksT* of the order of the De- 
bye energy, which is commonly used in literature.— We 
thus have 



scattering are practically the same for both samples. In 
fact, an equally good fit is obtained if the same T* is used 
for both curves. Note that the extracted T* ^ 200 K 
is indeed of the order of the Debye temperature 
of this material and thus leads to the conjecture that 
mostly acoustic phonons are involved in this scattering 
process.— Second, the spinon-defect scattering length Iq, 
which represents a lower bound for the low-T limit of 
^mag and which should significantly depend on the sam- 
ple's purity turns out to be drastically difi^erent for both 
cases. To be specific, we find Iq « 300 nm for the 2N com- 
pound and an extraordinary Iq ~ 1-6 iim for the 4N sam- 
ple, which correspond to more than 750 and 4100 lattice 
spacings, respectively. These findings provide a further 
confirmation of the above interpretation that, in both 
cases, Kinag is determined by the same spinon-phonon 
scattering process and that the difference between the 
two curves can be described by the different defect den- 
sity only. We mention that our results are consistent with 
recent data by T. Kawamata et al. — 

A major outcome of our study is the unambiguous 
identification of the extrinsic scattering processes as the 
only relevant ones. Intrinsic spinon-spinon scattering, on 
the other hand, plays no role in our analysis, even in the 
case of the very clean sample. The strong enhancement 
of Kmag upon reduction of the impurity amount thus ap- 
pears as the manifestation of ballistic heat transport of 
the underlying spin model, where Atmag is rendered finite 
by extrinsic scattering processes only. One might there- 
(1) fore speculate that Kmag of this material can be driven to 
much higher values in a perfect crystal. 

We point out that our analysis relies on a very sim- 
ple theoretical approach which was also successfully 
used in many other low-dimensional 5 = 1/2 spin sys- 
tems )2r2&i22ii2S which is surprising in view of the strong 
quantum nature of such systems. More sophisticated ap- 
proaches might lead to a deeper understanding of the 
magnetic heat transport in this system on a microscopic 
level. In this regard it is interesting to note that, in clean 
samples (i.e. with large Iq), = {AsT)~^ exp{T*/T) 



'mag ^ '0 



exp(T:/r) 



(2) 



which can be used to fit the data with Iq, Ag and T* 
(As describes the coupling strength) as free parameters. 
We find an excellent agreement between such fits and 
the experimental /mag, see Fig. [3l Inspection of the fit 
parameters^ yields two remarkable aspects which cor- 
roborate our previous qualitative findings. First, the pa- 
rameters As and T* which determine the spinon-phonon 



leads to k.^ w Kmag oc exp(T*/r) with T* - 200 K 
at high T, in agreement with the theory proposed 
by Shimshoni et ali^ However, we do not observe 
f^a = Kph 6xp(2T*/T) as expected in the same model. 
It seems worthwhile mentioning in this regard that the 
only slight enhancement of Ka observed upon increasing 
purity is quite unexpected. One might speculate that this 
is an indication of phonon scattering off the spin chains 
which in principle should be relevant i^i^^ 

To sum up, we have investigated the spinon heat con- 
ductivity Kmag of the antiferromagnetic S=l/2 Heisen- 
berg chain cuprate SrCu02 for standard (99%) and high 
(99.99%) purity. The higher purity leads to a drastic en- 
hancement of Kmag at low T and we find the up-to-present 
by far highest reported K^ag in the high-purity sample. 
For higher T, we provide clear-cut evidence that spinon- 
phonon scattering is the most relevant scattering which 
leads to a very efficient reduction of Kmag- An extreme 
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Figure 3: (color online) Magnetic mean free paths of SrCu02 
for different purities. The solid lines were calculated according 
to Eq.[2] The shaded area illustrates the uncertainty from the 
estimation of the phononic background. 



sensitivity of K,nag to impurities is present at low T , which 
implies that the spinon-defect scattering is dominating in 
this regime. A simple analysis reveals a remarkable lower 
bound for the low-temperature limit of the spinon mean 
free path /mag of more than a micrometer. Our results 
therefore suggest that Kmag is only limited by extrinsic 
scattering processes which appears as the manifestation 
of the ballistic nature of heat transport in the S = 1/2 
antiferromagnetic Heisenberg chain. 
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(1.56 ± 0.16) pm, T„*4jv = (204±11)K, Asan = 
(58.6 ±5.4) 10"^V/K and Zo,2iv = (305 ± 5) nm, 
TI2N = (217 ± 3) K, As,2N = (78 ±1) 10~"m/K, re- 
spectively. Setting T*2N ~ 7"u,4jv = 204 K after fit- 
ting the 4N data gives Zo,2JV = (320 ± 13) nm and 
As,2N = (72 ± 5) 10"^®m/K for 2N. The errors account for 
the accuracy of the fit. It is also possible to obtain a good 
fit with the same ^4^ for the 4N and 2N cases. However, in- 
dividual As account for errors in the absolute value, while 



fixing the energy scale by seems physically justified. 
Modeling Imag with l^p oc exp(r*/r) which accounts for 
the alternative scenario of spinous scattering off optical 
phonons22 results in a fit of similar quality with comparable 
Zo and T* ~ 300 K. 
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